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Abstract. A new picture of pulsar high-energy emission is proposed that is different 
from both the traditional polar cap and outer gap models, but combines elements of 
each. The slot gap model is based on electron acceleration along the edge of the open 
field region from the neutron star surface to near the light cylinder and thus could 
form a physical basis for the two-pole caustic model of (Dyks & Rudak, 2003). Along 
the last open field line, the pair formation front rises to very high altitude forming 
a slot gap, where the accelerating electric field is unscreened by pairs. The resulting 
radiation features both hollow cones from the lower-altitude pair cascades, seen at 
small viewing angles, as well as caustic emission on the trailing-edge field lines at 
high altitude, seen from both poles at large viewing angle. The combination of the 
small solid angle of slot gap emission (<C 1 sr) with a high probability of viewing the 
emission predicts that more gamma-ray pulsars could be detected at larger distances. 
In this picture, many of the positional coincidences of radio pulsars with unidentified 
EGRET sources become plausible as real associations, as the flux predicted by the 
slot gap model for many of the pulsars would provide the observed EGRET source 
flux. The expected probability of seeing radio-quiet gamma-ray pulsars in this model 
will also be discussed. 
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1. Introduction 

The number of rotation-powered pulsars with detected emission at X- 
ray and 7-ray energies has been steadily growing (Thompson, 2001; 
Becker & Aschenbach, 2002). Although the theory of pulsar acceleration 
and high-energy emission has been studied for over thirty years, the 
origin of the pulsed non-thermal emission is a question that remains 
unsettled. The observations to date have not been able to clearly dis- 
tinguish between an emission site at the magnetic poles (Daugherty & 
Harding, 1996) and emission in the outer magnetosphere (Cheng, Ho 
& Ruderman, 1986; Hirotani & Shibata, 2001). In the case of polar 
cap (PC) models, while the energetics and pair-cascade spectrum have 
had success in reproducing the observations, the predicted beam size 
of radiation emitted near the neutron star (NS) surface is too small to 
produce the wide pulse profiles that are observed (Thompson, 2001). 
However, Arons (1983) first noted the possibility of a high-altitude 
acceleration region or "slot gap" along the last open field line where 
the electric field is lower. The slot gap forms because the pair for- 




© 2008 Kluwer Academic Publishers. Printed in the Netherlands. 



SlotGap_HK.tex; 2/02/2008; 12:59; p.l 



2 



Harding & Muslimov 



mation front, above which the accelerating field is screened, occurs 
at increasingly higher altitude as the magnetic colatitude approaches 
the edge of the open field region (Arons & Scharleman, 1979 [AS79]). 
We have re-examined the slot-gap model (Muslimov & Harding, 2003 
[MH03]) with the inclusion of two new features: 1) the acceleration 
due to inertial-frame dragging (Muslimov & Tsygan, 1992) and 2) the 
additional decrease in the electric field near the boundary at the edge 
of the polar cap due to the narrowness of the slot gap. These two 
effects combine to enable acceleration to altitudes approaching the light 
cylinder in the slot gap at all azimuthal angles around the polar cap. 
These features result in the production of a larger high-energy emission 
beam with small solid angle, both favorable for producing high fluxes 
for 7-ray pulsars. 



2. Formation of the Slot Gap 

In the space-charge limited flow acceleration model AS79, electrons are 
freely emitted from the neutron star surface near the magnetic poles 
and accelerated in steady-state. When the electrons achieve a sufficient 
Lorentz factor they radiate curvature and inverse Compton photons 
that can produce electron-positron pairs in the strong magnetic field. 
Some of the positrons turn around and accelerate downward toward the 
neutron star, while the electrons accelerate upward. This pair polariza- 
tion screens the electric field above a pair formation front (PFF). These 
models assume a boundary condition that the accelerating electric field 
and potential are zero at the last open field line. Near the boundary, the 
electric field is lower and a larger distance is required for the electrons 
to accelerate to the Lorentz factor needed to radiate photons energetic 
enough to produce pairs. The PFF thus occurs at higher and higher 
altitudes as the boundary is approached and curves upward, approach- 
ing infinity and becoming asymptotically parallel to the last open field 
line. If the electric field is effectively screened above the PFF, then a 
narrow slot surrounded by two conducting walls is formed (see Figure 
1). Within this slot gap, the electric field is further screened by the 
presence of the second conducting boundary and acceleration occurs 
at a reduced rate. Pair production and pair cascades therefore do not 
take place near the neutron star surface in the slot gap, as do the pair 
cascades along field lines closer to the magnetic pole (core), but occur 
at much higher altitudes. 

There are several important differences between the revised slot-gap 
model of MH03 and the original slot-gap model of AS 79. The inclusion 
of general relativistic frame dragging enables particle acceleration on 
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Figure 1. Schematic illustration of polar cap geometry, showing the outer boundary 
of the open field line region (where Eu = 0) and the curved shape of the pair 
formation front (PFF) which asymptotically approaches the boundary at high al- 
titude. The slot gap exists between the pair plasma which results from the pair 
cascades above the PFF and the outer boundary. A narrow beam of high-energy 
emission originates from the low-altitude cascade on field lines interior to the slot 
gap. A broader, hollow-cone beam originates from the high-altitude cascade above 
the interior edge of the slot gap. A£ SG is the slot gap thickness (see text) and 9q,sg 
is the colatitude at the center of the slot gap. 



both "favorably" curved (curving toward the rotation axis) and "unfa- 
vorable" curved (curving away from the rotation axis) field lines and 
also at all inclination angles. MH03 also consider the radiation from 
pair cascades occurring along the interior edge of the slot gap. The 
cascade radiation emission beam from the slot gap is thus a full hollow 
cone centered on the magnetic axis. A narrower emission beam from 
field lines interior to the slot gap will form a core component of pairs 
and high-energy emission. 

Since the potential in the slot gap is uncreened, electrons on field 
lines which thread the slot gap will continue accelerating to very high al- 
titudes. Muslimov k, Harding(2004, [MH04]) have derived the potential 
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and accelerating electric field at high altitude in the extended slot gap 
by considering the effect of deviation of stream lines from the magnetic 
field lines of a static dipole. In the absence of this effect, the difference 
between the actual charge density along field lines and the corota- 
tion (or Goldreich-Julian) charge density, which grows with altitude 
above the polar caps, would become comparable to the Goldreich-Julian 
charge density itself, a situation which cannot be supported since it 
would disrupt the steady-state charge flow within the magnetic flux 
tube. However, long before this would occur, the drift of electrons in 
the slot gap across the magnetic field will largely damp the growth of 
the charge deficit and the large electric field which would be induced 
perdendicular to the magnetic field. The residual parallel electric field 
is small and constant, but still large enough at all altitudes up to nearly 
the light cylinder to maintain a flux of electrons with Lorentz factors 
exceeding 10 7 in the slot gap. 

MH04 matched the high-altitude slot gap solution for the parallel 
electric field to the solution found at lower altitudes (MH03). The re- 
sult, for most inclination angles, is a continuously accelerating field from 
the neutron star surface to near the light cylinder along the last open 
field lines. The derived properties of the extended slot gap closely match 
the geometrical properties required for the two-pole caustic model of 
Dyks & Rudak (2003). Such a model can reproduce the double-peaked 
profiles seen in many 7-ray pulsars like the Crab, Vela and Geminga. 
MH04 demonstrated through numerical simulation of slot gap acceler- 
ation and radiation, that the extended slot gap radiation will produce 
the caustic peaks and pulse profiles similar to those of the two-pole 
caustic model. 

The existence of a slot gap requires a dense enough pair plasma on 
interior field lines to fully screen the parallel electric field, so that the 
inside walls of the gap have vanishing electric field. Harding & Mus- 
limov(2002, HM02) found that pairs produced by CR are always able 
to fully screen the E\\ whereas pairs from ICS were only able to screen 
locally, if at all. Therefore we can conclude that only those pulsars 
which can produce CR pairs will have slot gaps. Figure 2 reproduces the 
death lines for CR pairs in P-P space from HM02, showing the region of 
pulsar parameter space required for slot gap formation. Generally, the 
younger pulsars, with ages less than 10 7 yr and with higher magnetic 
fields fall into this region, although one or two millisecond pulsars may 
also have slot gaps. 
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Figure 2. Radio pulsars from the ATNF catalog plotted on the P-P diagram. The 
gray dots are radio pulsars, solid squares are EGRET pulsars, open triangles are 
radio pulsars in EGRET sources error circles. The gray line is the curvature radiation 
pair death line. 

3. Slot Gap Energetics 

The electrodynamics of the slot gap is primarily dependent on a single 
parameter, the slot gap width, A£ SG . The ratio of the electric field 
in the slot gap to the electric field in the core region of the PC is 
E SG I .E^orc oc A£;? G / 4 and the luminosity of particles accelerated in the 
slot gap is .Lprim oc A^ g Lsd where Lsd is the spin-down luminosity 
of the pulsar (for full details see Muslimov & Harding, 2003). One can 
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Figure 3. Observed flux above 1 keV, <3? 7 , times distance squared (from Thompson 
(2001)) (solid circles) and theoretical values of specific high-energy luminosity from 
the slot gap, e 7 -I/ pr im/£2 7 from eq. (4) (upside-down triangles) vs. spin-down lumi- 
nosity for known 7-ray pulsars. An efficiency of e 7 = 0.3 was assumed. Also A = 0.1, 
»7 7 = 3 and the stellar parameters Rs = 1.6 and Lts = 4 were used. 



estimate the width of the slot gap as the magnetic colatitude where the 
variation in height of the curvature radiation PFF becomes comparable 
to a fraction A of the stellar radius R, or 

ro.2ifa(AB 13 )- 4 /fj-3A B 12 <4A 
§SG ~ 1 0.1 PosiXB^r 4 / 7 ^ 7 Bu > 4.4 

I45 = //10 45 g-cm 2 , I is the neutron star moment of inertia, Pn.l = 
P/0.1 s and B12 = Bq/W 12 G are the neutron star rotation period and 
surface magnetic field. Here, A£ SG is in units of the polar cap half-angle, 
#0 = sin _1 (27ri?/Pc) 1//2 . The emission solid angle of radiation from the 
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slot gap can be estimated by integrating over the thin annulus defined 
by the slot gap width (Eqn [1]). 

9 

Msg ~ -7r6» ?7A£ 5G ster, (2) 

where rj = r/R is the dimensionless radius of emission. Electrons accel- 
erating in the slot gap will radiate curvature-radiation 7-rays, becoming 
radiation-reaction limited at Lorentz factors 



7 « 3 • 10 7 



K0.15B12— ^ SG |cosx| 

Vic 



1/4 



(3) 



where k 0A5 = k/0.15 » I45/RI, u sg = °- 25A ^L' ^6 = R/W 6 cm 
is the neutron star radius, \ is the pulsar inclination angle and r\[ c is 
the dimensionless light cylinder radius. Based on the luminosity of the 
primary electrons and the above solid angle estimate, we may derive 
the quantity, 

^SG 

= 3 x 10 34 e 7 L%,P^R l J'\- l \-*lX! n ^0.15 cos 2 X ergs" 1 ster^) 

where L pr im is the luminosity in primary electrons accelerated in the 
slot gap, e 7 is the radiation efficiency and L s d,35 = -^sd/10 35 ergs -1 is 
the spin-down luminosity. The above expression for L 7 (f2 7 ) is equiv- 
alent to the observed quantity <I> 7 d 2 , where <3? 7 is the high-energy 
bolometric flux observed at the Earth, and d is the distance to the 
pulsar. Figure 3 shows the observed (solid circles with error bars) and 
theoretical (upside-down triangles) values of <E> 7 (i 2 as a function of spin- 
down luminosity, L s d for the 10 known 7-ray pulsars. The theoretical 
values are calculated for the parameters e 7 = 0.3 and A = 0.1 (see 
eq. [4]). Note that parameter e 7 can range from 0.2 to 0.5 in cascade 
simulations, and 77 = 3. In Figure 3 the dashed line represents the limit, 
where the spin-down luminosity is radiated into the unit solid angle, 
i.e. where $ 7 <i 2 = L s d/1 ster. One can see that there is good agreement 
for most high-energy pulsars except several of the pulsars, Geminga 
and PSR B0656+14, having low L sd , and for J0218+4232, which is 
a millisecond pulsar. These pulsars are near or below the curvature 
radiation pair death line (see Figure 2 and HM02), and therefore have 
either very wide slot gaps or no slot gaps at all. All other high-energy 
pulsars depicted in Figure 3 are well above the curvature radiation 
death lines and are expected to have slot gaps. 
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Figure 4- Comparison of fluxes of EGRET sources (white bars) and of associated 
radio pulsars (gray bars) as predicted by the slot gap model (Eqn [4] and Muslimov 
& Harding 2003) 



4. Slot Gap Emission and EGRET Unidentified Sources 



There are presently 53 radio pulsars in the error circles of 28 EGRET 
unidentified sources (Kramer et al., 2003; Grenier, 2004). Nearly all 
of these pulsars have been discovered by surveys such as the Parkes 
Multibeam (Manchester et al., 2001) and by deep pointed searched 
by Arecibo and GBT radio telescopes (Camilo, 2004) after the end of 
EGRET operation. The pulsars coincident with EGRET sources are 
plotted in P-P space in Figure 2. Most of them are young (r < 10 7 yr), 
high-field pulsars above the CR pair death line and many are expected 
to have slot gaps. There are unfortunately too few photons in any of the 
sources to do credible pulsation searches of the EGRET archival data, 
so such searches must await AGILE or GLAST. Many of the coincident 
EGRET sources contain multiple radio pulsars within the error circles. 
These are shown in Figure 4 with their observed average fluxes from 
the 3rd EGRET catalog (Hartman, 1999). 
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We have computed the 7-ray fluxes predicted by the slot gap model 
(Eqn [4]) for the EGRET source-coincident radio pulsars and the results 
are plotted in Figure 4. The fluxes of the pulsars, LsG(^)/d 2 , where d 
is the distance, are shown as dark bars alongside the coincident EGRET 
source fluxes shown as light bars. In about 18-22 out of the 28 sources, 
the predicted slot gap flux from either an individual pulsar or several 
pulsars combined could account for the EGRET source flux. Thus more 
than two thirds of these associations are physically plausible, making 
the pulsars viable counterparts for the EGRET sources. In contrast, 
the predicted fluxes from the standard polar cap model (e.g. HM02), 
assuming a solid angle of 1 sr., would be comparable to the EGRET 
source fluxes in only about 5 of the cases. 



5. Conclusions 



In the slot gap model, pulsar high-energy emission comes from high 
altitudes, between 2-3 stellar radii above the neutron star surface to 
nearly the light cylinder. Two resulting characteristics of radiation from 
the slot gap: small solid angle and a wide emission beam, combine to 
provide a significantly larger flux for 7-ray pulsars than in the stan- 
dard polar cap model. The small solid angle allows the radiation from 
7-ray pulsars to be seen by any given detector at larger distances, 
and many of the spatially coincident pulsars are at large distances, as 
determined from their dispersion measure. Thus, a larger number of 
the recently detected radio pulsars that are in or near EGRET source 
error circles become plausible candidates for the 7-ray sources. It is 
possible that half of the non-variable unidentified EGRET sources in 
the galactic plane are radio loud 7-ray pulsars. A population synthesis 
study of radio-loud and radio-quiet pulsars in the galaxy (Gonthier 
et al., 2004a; Gonthier et al., 2004b) have found that in the slot gap 
model one expects that many more radio-loud than radio-quiet pulsars 
are counterparts to EGRET sources. This result is due to the fact that 
the radio emission occurs along the same open field lines as the 7-ray 
emission, in contrast to outer gap models where the radio emission is 
assumed to originate from open field lines on the opposite pole of the 
neutron star from the high-energy emission. The ratio of radio-loud 
to radio-quiet 7-ray pulsars detected by EGRET, and eventually by 
GLAST, will be a good discrimator between emission models. 
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